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Abstract 

The Misner and Sharp approach to the study of gravitational col- 
lapse is extended to the dissipative case in, both, the streaming out 
and the diffusion approximations. The role of different terms in the 
dynamical equation are analyzed in detail. The dynamical equation is 
then coupled to a causal transport equation in the context of Israel- 
Stewart theory. The decreasing of the inertial mass density of the 
fluid, by a factor which depends on its internal thermodynamics state, 
is reobtained, at any time scale. In accordance with the equivalence 
principle, the same decreasing factor is obtained for the gravitational 
force term. Prospective applications of this result to some astrophys- 
ical scenarios are discussed. 
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1 Introduction 



Some years ago, Misner and Sharp [1] and Misner [2] provided a full account 
of the dynamical equations governing the adiabatic [1], and the dissipative 
relativistic collapse in the streaming out approximation [2]. 

The relevance of dissipative processes in the study of gravitational col- 
lapse cannot be over emphasized. Indeed, dissipation due to the emission 
of massless particles (photons and/or neutrinos) is a characteristic process 
in the evolution of massive stars. In fact, it seems that the only plausible 
mechanism to carry away the bulk of the binding energy of the collapsing 
star, leading to a neutron star or black hole, is neutrino emission [3]. 

In the diffusion approximation, it is assumed that the energy flux of 
radiation, as that of thermal conduction, is proportional to the gradient of 
temperature. This assumption is in general very sensible, since the mean free 
path of particles responsible for the propagation of energy in stellar interiors 
is in general very small as compared with the typical length of the object. 
Thus, for a main sequence star as the sun, the mean free path of photons 
at the centre, is of the order of 2 cm. Also, the mean free path of trapped 
neutrinos in compact cores of densities about 10^^ g. cm~^ becomes smaller 
than the size of the stellar core [4, 5]. 

Furthermore, the observational data collected from supernova 1987A in- 
dicates that the regime of radiation transport prevailing during the emission 
process, is closer to the diffusion approximation than to the streaming out 
limit [6]. 

However in many other circumstances, the mean free path of particles 
transporting energy may be large enough as to justify the free streaming 
approximation. Therefore it is advisable to include simultaneously both lim- 
iting cases of radiative transport, diffusion and streaming out, allowing for 
describing a wide range of situations. 

In a recent work [7] we have studied the effects of dissipation, in both 
limiting cases of radiative transport, within the context of quasi-static ap- 
proximation. This assumption is very sensible because the hydrostatic time 
scale is very small for many phases of the hfe of a star. It is of the order 
of 27 minutes for the sun, 4.5 seconds for a white dwarf and 10"'' seconds 
for a neutron star of one solar mass and 10km radius [8]. However, during 
their evolution, self-gravitating objects may pass through phases of intense 
dynamical activity, with time scales of the order of magnitude of (or even 
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smaller than) the hydrostatic time scale, and for which the quasi-static ap- 
proximation is clearly not reliable, e.g., the collapse of very massive stars [9], 
and the quick collapse phase preceding neutron star formation, see for exam- 
ple [10] and references therein. In these cases it is mandatory to take into 
account terms which describe departure from equilibrium, i.e. a full dynamic 
description has to be used. 

Thus our primary task here consists in extending the Misner dynamical 
equations as to include dissipation in the form of a radial heat flow (besides 
pure radiation). This is presented in Section 3. Then in the following Sec- 
tion, the resulting dynamical equation is coupled to the transport equation 
obtained in the context of the Miiller-lsrael-Stewart theory [11, 12]. 

After doing that we show that the effective inertial mass density of a fluid 
element reduces by a factor which depends on dissipative variables. This 
result was already known (see [13] and references therein), but considered to 
be valid only, just after leaving the equilibrium, on a time scale of the order 
of relaxation time. The novelty here is, on the one hand, that such reduction 
of the effective inertial mass density is shown to be valid at an arbitrary time 
scale, and on the other, that the "gravitational force" term in the dynamical 
equation is also reduced by the same factor, as expected from the equivalence 
principle. Prospective applications of this result to astrophysical scenarios 
are discussed at the end. 

In the next section the fleld equations, the conventions, and other useful 
formulae are introduced. 

2 The energy— momentum tensor and the field 
equations 

In this section we provide a full description of the matter distribution, the 
line element, both, inside and outside of the fluid boundary and the fleld 
equations this line element must satisfy. Since we are going to follow closely 
the Misner approach [2] we shall use comoving coordinates (for a description 
of gravitational collapse in non-comoving coordinates, see [14] and references 
therein) . 
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2.1 The interior spacetime 

We consider a spherically symmetric distribution of collapsing fluid, which, 
for sake of completeness, we assume to be locally anisotropic, undergoing 
dissipation in the form of heat flow and free streaming radiation, bounded 
by a spherical surface E. For such system the energy-momentum tensor is 
given by 

T^P = (/i + P^'^V^ + P^g''^ + {Pr - Pi)x'"X^ + (fy^^ + ^"(f + e/"/^, (1) 

where, ;U is the energy density, Pr the radial pressure, Pj_ is the tangential 
pressure, e is the radiation density, is the four velocity of the fluid, (f is 
the heat flux, x° is a unit four vector along the radial direction and P is a 
null four vector. These quantities have to satisfy 

V"K. = -1, V"g, = 0, x"Xa = l, X"K. = 0, /X = 0. (2) 

We assume the interior metric to S to be comoving, shear free for simplicity, 
and spherically symmetric, accordingly it may be written as 

and hence 

V'^^A-^b'^, = r ^ A-^8'^ + B-^8'i, = ^"^<^f , (4) 

where g is a function of t and r and we have numbered the coordinates — t, 
= r, x'^ = 9 and x^ — (p. Now the Einstein's fleld equations become with 
the help of (1-4) 
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where the dot and prime stand for differentiation with respect to t and r. 
The rate of expansion © = V'^-a of the fluid sphere is given, from (3) and 
(4), by 

© = 34., (9) 



and from (6) we have 



AB' 



8n{qB + e)B = -9'. 



(10) 



If g > and e > 0, then from (10) we have 9' > meaning that, if the 
system is collapsing 9 < 0, g and/or e decrease the rate of collapse towards 
the outer layers of matter. If g = and e = from (10) 9' = 0, which means 
that the collapse is homogeneous. 

The mass function m(t, r) of Cahill and McVittie [15] is obtained from 
the Riemann tensor component i?23^^ and is for metric (3) 



m(i, r) 



7, -"-23 



^3 ^^1 



12 



r^B 



- r^B'. 



(11) 



2.2 The exterior spacetime 



The exterior spacetime to E of the collapsing body is described by the out- 
going Vaidya spacetime which models a radiating star and has metric 



2m{v) 



dv^ - 2dvdp + p^de"^ + sin^ ed(j)'^), 



(12) 



where m, the total mass inside S, is a function of the retarded time v. The 
surface E described by the comoving coordinate system (3) is r = = 
constant, while in the non comoving coordinate system (12) is p = Pt,{v). 
Matching the interior spacetime (3) with source (1) to the exterior spacetime 



5 



(12) by using Darmois junction conditions we obtain 



(13) 



( 



2 ^2 2 B 



3 BB^ r3 



{qB + e)s 




(14) 



(15) 



(16) 




) 



1/2 



(17) 



s 



where L is defined as the total luminosity of the collapsing sphere as measured 
on its surface and is given by 



and where 

is the total luminosity measured by an observer at rest at infinity. The result 
(13) represents the continuity of the radial flux of momentum across E in 
which only the heat flow q appears. However, for the total radiation leaving 
S (14) the radiation e contributes as well as q. Although it might seem to 
be obvious, it is perhaps important to stress that the radiation e has the 
same null property associated to the exterior null radiation that it produces, 
while the heat flux q, producing the exterior null radiation too, is not a 
null flux. Relation (15) is the equality of the proper radii as measured from 
the perimeter of the spherical surface S in both frames (3) and (12). The 
expression for the total mass (16) is the corresponding mass function [15] 
given by (8). The relationship of proper times measured on E with both 
frames (3) and (17) is given by (17). 




(18) 
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3 Dynamical equations 



For studying the dynamical properties of the field equations and following 
Misner and Sharp, let us introduce the proper time derivative Dt given by 



1^ 



Dt = 4^- (20) 
Then using (20) we can describe the velocity U of the collapsing fiuid as 



U — rDfB < {in the case of collapse), 
then (11) can be rewritten as 



(rB)' 
B 



1 + - 



2m{t, r) 
rB 



1/2 



E. 



(21) 



(22) 



The right hand side of (22) can be interpreted as being the energy density E 
of a collapsing fiuid element. Next, by taking the proper time derivative of 
(11) we obtain 



Dtm 



,BBB 

\ 
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B 



.BAB^ 



— r 



A^ 



r^BB'^ .B'B' oB' 

H r r — . 

' AB A 



2 AB^ 

Considering (6) and (7) we can rewrite (23) as 

.B^B 



(23) 



Dtui — —Att 



[Pr + e)r' 



A 



+ {qB + e)r^B{rB) 



and with (21) and (22) it becomes 

Dtm = 47r [-(P, + e)U - {qB + e)E] {rBf, 



(24) 



(25) 



which gives the rate of variation of the total energy inside a surface of radius 
rB. In the right hand side of (25) {Pr + ^)\U\ (in the case of collapse f/ < 0) 
increases the energy inside rB through the rate of work being done by P^ 
and the induction field produced by e and already observed by Misner [2]. 
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Clearly here the heat flux q does not appear since it does not produce an 
induction field. The second term —{qB + e)E is the matter energy leaving 
the spherical surface. 

Another proper derivative that helps us to study the dynamics of the 
collapsing system is the proper radial derivative Dji, where 



R^rB, 



(26) 



constructed from the radius of a spherical surface, as measured from its 
perimeter inside E, being 

1 d 



R'dr' 



(27) 



Then by taking the proper radial derivative (27) of (11) we obtain 
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Considering (5) and (6) then (28) becomes 
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(28) 



Djim — Att 



li + e+{qB + e) 



rBB 
{rByA 



{rBf 



and with (21) and (22) we finally have 



DRim = 4:7T 



li + e + {qB + e) 



{rBf 



(29) 



(30) 



This expression gives the total energy entrapped between two neighboring 
spherical surfaces with respect to proper radius inside the fiuid distribution. 
The first term on the right hand side of (30) ;U + e is due to the energy 
density plus the induction field and no heat fiux appears. The second term 
{qB + e)U/E is negative (in the case of collapse) and measures the out flux 
of heat and radiation. 

Finally, we can obtain the acceleration DfU of a collapsing particle inside 
E. In order to do that we start from (7) and (11) which allows us to write 

A' (rB)' 



m + A7r(Pr + e){rBy 



(rB)-' + 



A B^ 



(31) 
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Calculating the r component of the Bianchi identities, tI^.^ — 0, from (1) 
wc obtain 




Substituting the expression A' /A from (32) into (31) and considering (11), 
(20), (26) and (27) we obtain 



The first term with square brackets in the right hand side of (33) repre- 
sents the gravitational force. It shows that the gravitational force acting on a 
particle has a Newtonian part with m and a purely relativistic gravitational 
contribution due to and e. The second term in square brackets represent 
the hydrodynamical forces. It consists of the usual pressure gradient term 
(including the contribution of the radiation to the pressure) Dfi{Pr + e) < 
counteracting collapse, and the anisotropic force term e + Pr — P± which can 
be positive or negative thus, respectively, accelerating more or counteracting 
the rate of collapse. In these two terms the appearance of e is due to the 
contribution of radiation to the total energy density and radial pressure. The 
last term in square brackets contains the specific contribution of dissipation 
to the dynamics of the system. The first term within this bracket is pos- 
itive {U < 0) showing that the out flux oi q > and e > diminish the 
total energy inside the collapsing sphere thereby reducing the rate of col- 
lapse. It is interesting to observe the different effects that q and e have on 
the dynamical behaviour of the collapsing fluid. The heat flux q helps only 
to slow down the rate of collapse by diminishing the energy inside the fluid 




-E {5qB + Ae)- + BDtq + Dte . 



(33) 



Equation (33) has the "Newtonian" form 



Force — Mass density x Acceleration 



(34) 
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sphere by producing an exterior outflowing radiation. On the other hand, 
the radiation density e behaves not only in a similar way as q by diminishing 
the energy of the collapsing sphere through the exterior outflow of radiation, 
but contributes too as an induction field to the gravitational energy, first 
observed by Misner and collaborators [2, 16], and contributes to the radial 
pressure P^. The effects of Dfe have been discussed in detail in [2]. Thus it 
remains to analyse the effects of D^q, this will be done in the next section 
after introducing the transport equation. 

Before coming to the next section, we observe that from (33) the limit 
of hydrostatic equilibrium when U = 0, g = and e = can be achieved, 
producing 

which is just the generalization of the TOV equation for anisotropic fluids 
[17], obtained in comoving coordinates in [18] while studying dynamical in- 
stability for radiating anisotropic collapse. 



4 Transport equation and its consequences 

As we mentioned before we shall use a transport equation derived from the 
Miiller-Israel-Stewart second order phenomenological theory for dissipative 
fluids [11, 12]. 

Indeed, it is well known that the Maxwell-Fourier law for the radia- 
tion flux leads to a parabolic equation (diffusion equation) which predicts 
propagation of perturbation with infinite speed (see [19]-[22] and references 
therein). This simple fact is at the origin of the pathologies [23] found in 
the approaches of Eckart [24] and Landau [25] for relativistic dissipative pro- 
cesses. 

To overcome such difficulties, different relativistic theories with non-vanishing 
relaxation times have been proposed in the past [11, 12, 26, 27]. The impor- 
tant point is that all these theories provide a heat transport equation which 
is not of Maxwell- Fourier type but of Cattaneo type [28] , leading thereby to 
a hyperbolic equation for the propagation of thermal perturbation. Thus the 
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corresponding transport equation for the heat flux reads [21] 

Th''''V^qfS.,, + q'' = -Kh'^^{Tp + Ta0)--KT'i^--^j q\ (36) 

where h^^ is the projector onto the three space orthogonal to V^, n denotes 
the thermal conductivity, and T and r denote temperature and relaxation 
time respectively. Observe that, due to the symmetry of the problem, equa- 
tion (36) only has one independent component, which may be writtten as: 

T{qB)B + qAB^ = -k{TA)' 1— j (37) 

Which, using (20) and (31), becomes 
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(38) 



We can couple the transport equation in the form above (38) to the 
dynamical equation (33), in order to bring out the effects of dissipation (in 
the diffusion approximation) on the dynamics of the collapsing sphere. With 
that purpose, let us replace (38) into (33) (putting e = 0), then we obtain 
after some rearrangements 



ill + Pr){l - a)DtU = Fgrav{l - «) + Fhyd + 



EkT' EqB KET^qB ( r 
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2At 
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(39) 



Where Fg^av and F^yd are defined by 



Fgrav = -(a* + -Pr) m + AttP^R^ 



1 



(40) 



and 



hyd 



-E' 



DnPr + 2- 



P.-Pl 

R 



(41) 
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where a is given by 

« = -r^^- (42) 

We can now analyze the overall effects of dissipation (in the diffussion 
approximation) on the evolution of the collapsing sphere. 

First of all observe that as a tends to 1, the effective inertial mass density 
of the fluid clement tends to zero. This effect was known (sec [13] and 
references therein), but only to be valid just after the system abandons the 
equilibrium, on a time scale of the order of (or smaller than) the relaxation 
time. Here we see that it is present all along the evolution. Furthermore 
we see that Fgrav is also multiplied by the factor (1 — a). Indicating that 
the effective gravitational attraction on any fluid element decreases by the 
same factor as the effective inertial mass (density). Which of course is to be 
expected, from the equivalence principle. It is also worth mentioning that 
F^yd is in principle independent (at least explicitly) on this factor. 

Next observe that the third and the fourth terms, on the right hand side 
of (39), arc of oposite sign and of the same order of magnitude (at least in 
the case of not too strong gravitational field). Finally, the sign and the order 
of magnitude of the fifth term is clearly, model dependent. Furthermore this 
term is absent in the "truncated" version of the theory (see [29]). 

With these comments above in mind, let us imagine the following situ- 
ation: A collapsing sphere evolves in such a way that the value of a keeps 
increasing and approaches the critical value of 1. As this process takes place, 
the ensuing decreasing of the gravitational force term would eventually lead 
to a change of the sign of the right hand side of (39). Since that would hap- 
pen for small values of the effective inertial mass density, that would imply 
a strong bouncing of the sphere, even for a small absolute value of the right 
hand side of (39). At the origin of this effect, of course, is the equivalence 
principle, according to which the inertial mass and the passive gravitational 
mass are equal. It is also worth noticing that changes of inertial mass due to 
different physical phenomena, are familiar in other branches of physics. Thus 
for example the effective inertial mass of an electron moving under a given 
force through a crystal, differs from the value corresponding to an electron 
moving under the same force in free space, and may even become negative 
(see [30]). 

For this picture to be physcally meaningful one should first answer to the 
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following questions: 



• How close may a approach the critical value? 



• In what physical scenarios one could expect values of a close to the 
critical value? 

Since these questions are related to each other, we shall answer to them 
simultaneously. 

First of all it should be mentioned that from the analysis of stability 
and causality in dissipative relativistic fluids [23], it follows that causality 
and hyperbolicity (which imply stability) require for dissipative viscous free 
systems 



r > 



+ 



n c„ 



r > 



kT 




and 



T > 



K 



2apT 



(43) 



(44) 



(45) 



_Kt(/U +p) 

where n and Cg denote the particle density and the sound speed, Cp and are 
the specific heat at constant pressure and volume, kt is the thermal expansion 
coefficient and ctp the isothermal compressibility. These expressions are found 
from equations (146-148) in [23], taking the limit /3o, 00 and ccj = 

(this method was applied in [21] to the case in which only bulk viscous 
perturbations were present). It should be kept in mind that the conditions 
above, are obtained within a linear perturbative scheme. 

Obviously, condition (43) is violated at the critical point (in fact it is 
violated, slightly below it). However as we shall see below, it is not difficult 
to find physical conditions for which the numerical values of variables entering 
in the definition of a lead to a = 1. Therefore, the relevant question is: Can 
a physical system actually reach the critical point? If the answer to this 
question is negative, then it should be explained how a given system avoids 
the critical point. Since, as mentioned before, numerical values of k, T, t, 
/J,, and Pr, leading to a ~ 1 may correspond to a non very exotic scenario. 
On the other hand, a positive answer seems to be prohibited by causality 
and stability conditions. However, we shall conjecture here that this might 
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not be the case. In fact, the vanishing of the effective inertial mass density 
at the critical point, indicates that linear approximation is not valid at that 
point. So it seems that the behaviour of the system close to the critical point 
cannot be studied with a linear perturbative scheme, in which case it might 
be possible for a given system to attain the critical point. Furthermore, in the 
general case (including viscosity) it may happen that causality breaks down 
beyond the critical point [31]. Thus, it appears that there exist situations 
where a given physical system may attain the critical point and even go 
beyond it. 

Indeed, condition a ^ 1 can be accomplished in non very exotic systems. 
One of them is an interacting mixture of matter and neutrinos, which is a 
well-known scenario during the formation of a neutron star in a supernova 
explosion. In this case the heat conductivity coefficient is given by [32, 33] 

K = ^bT\ (46) 

where r is the mean collision time and b = 7iVjya/8, with Ni, the number 
in neutrino flavors and a the radiation constant. Assuming that the two 
viscosity coefficients vanish, and since p should be not larger than /i, then 

-~ — . (47) 

Using usual units, the critical point is overtaken if 

where we have adopted r ~ t«;, A^,^ = 3, T is in Kelvin and /j, is given in g 
cm~^. The values of temperature, for which a ~ 1, are similar to the expected 
temperature that can be reached during hot collapse in a supernova explosion 
[34, § 18.6] 

Also, we may speculate that a may increase substantially (for non-negligible 
values of r) in a pre-supernovae event 

Indeed, at the last stages of massive star evolution, the decreasing of 
the opacity of the fluid, from very high values preventing the propagation of 
photons and neutrinos (trapping [4]), to smaller values, gives rise to radiative 
heat conduction. Under these conditions both k. and T could be sufficiently 
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large as to imply a substantial increase of a. Indeed, the values suggested in 
[35] {\k\ « 10^^; [T] ^ 10^3; [r] ^ 10"^; [p] ^ 10^^ in c.g.s units and Kelvin 
) lead OL'^l. The obvious consequence of which would be to enhance the 
efficiency of whatever expansion mechanism, of the central core, at place, 
because of the decreasing of the inertial mass density and the gravitational 
force term. At this point it is worth noticing that the relevance of relaxational 
effects on gravitational collapse has been recently exhibited and stressed (see 
[36] , and references therein) 

5 Conclusions 

Following the scheme developped by Misner and Sharp, we have established 
the set of dynamical equations governing the evolution of collapsing dissipa- 
tive spheres, taking into account, both, the free-streaming and the difussion 
approximation. We have further coupled the dynamical equation with a heat 
transport equation obtained from the Miiller-Israel- Stewart theory. The re- 
sulting equation brings out the relevance of a critical point (a = 1) for 
which, both the effective inertial mass density and the gravitational frorce 
term vanish. We have shown that in principle that critical point may be at- 
tained under acceptable physical conditions (e.g. a supernova scenario) and 
have speculated about the possibility that in that collapsing sphere 

bounces. Of course the eventual application of this model to a supernova, 
would require much more details about the astrophysical settings. 
Before concluding we would like to make the following remarks: 

1. It should be noticed that the appearance of the factor 1 — a in the 
inertial mass density and the gravitational force term, is related to the 
first term on the left of equation (36). But this is the term which 
introduces causality in the transport equation and therefore, is to be 
expected in any causal theory of dissipation. Accordingly our main 
result is also expected to hold for a general family of theories which 
includes the Miiller-Israel-Stewart theory. However, whereas the mere 
appearance of the factor 1 — a in the dynamical equation is a very 
general result, the possibility of reaching the critical point and the 
physical consequences derived from that, will depend on the specific 
theory of dissipation to be adopted. In this same line of arguments. 
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it should be noticed that the transport equations in the Miiller-Israel- 
Stewart theory coincide with those obtained from the kinetic theory. 
Accordingly the critical point should also appear in a kinetic theory 
approach. 

2. Observe the formal similarity between the critical point and the equa- 
tion of state for an inflationary scenario = —Pr) without dissipation. 
This kind of equation of state has been recently proposed to describe 
the interior of a cold compact object without event horizons, and which 
would represent an alternative to black holes [37] . One could speculate 
with the possibihty that such interior could be described instead, by a 
dissipative fluid with a tending to 1. 

3. It should be clear that the analysis presented here depends strictly on 
the validity of the diffusion approximation, which in turn depends on 
the assumption of local thermodynamical equilibrium (LTE). There- 
fore, only small deviations from LTE can be considered in the context 
of this work. Thus when we state that the effective incrtial mass den- 
sity decreases by the factor (1 — a), at all time scales, this means at 
time scales within which the system is not very far from LTE. 

4. For the sake of completeness we have considered an anisotropic fluid (in- 
stead of an isotropic one, Pr = -Pi), leaving the origin of such anisotropy 
completely unspecifled. As it is apparent, anisotropy does not affect 
the most important result obtained here (e.g. the existence of the criti- 
cal point). However, should anisotropy be related to viscosity, then for 
consistency the anisotropic pressure tensor should be subjected to the 
Israel-Stewart causal evolution equation for shear viscosity. 
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